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Abstract 

A quantal guiding center theory allowing to systematically study the separa- 
tion of the different time scale behaviours of a quantum charged spinning particle 
moving in an external inhomogeneous magnetic filed is presented. A suitable set 
of operators adapting to the canonical structure of the problem and generalizing 
the kinematical momenta and guiding center operators of a particle coupled to a 
homogenous magnetic filed is constructed. The Pauli Hamiltonian rewrites in this 
way as a power series in the magnetic length Ib = yUc/eB making the problem 
amenable to a perturbative analysis. The first two terms of the series are explic- 
itly constructed. The effective adiabatic dynamics turns to be in coupling with a 
gauge filed and a scalar potential. The mechanism producing such magnetic-induced 
geometric-magnetism is investigated in some detail. 



short title: A Quantal Guiding Center Theory 
PACS: 03.65.-W, 02.40+m, 52.20. Dq, 02.90.+p 



1 Introduction 



The motion of a charged particle in a strong inhomogeneous magnetic filed is a nontrivial 
problem displaying a variety of very interesting phenomena ranging from chaos to phase 
anholonomy. Being of utmost importance in plasma physics, expecially in the study of 
magnetic confinement, the subject has been worked out in great detail in classical mechan- 
ics with special attention to phenomenological implications as well as to formal aspects. 
The canonical structure of the problem, in particular, has been deeply investigated only 
in a relatively recent time by R.G. Littlejohn [|I| revealing the appearance of geometry 
induced gauge structures in the adiabatic motion of classical charged particles. Very few, 
on the other hand, is known about the behaviour of quantum particles in strong inho- 
mogeneous magnetic fields, the reason being essentially that the techniques developed 
for classical mechanics do not easily generalize to the quantum context. Some work has 
been done for neutral spinning particles by M.V. Berry p[, Y. Aharonov & A. Stern 0] 
and R.G. Littlejohn Sz S. Weigert in connections with geometrical phases, whereas a 
quantum treatment for charged spinning particles is still missing. It is the purpose of 
this paper to present what may be probably called a quantal guiding center theory in 
which the coupling between the spin and spatial degrees of freedom of a quantum charged 
spinning particle moving in a strong inhomogeneous magnetic field is systematically taken 
into account. This allows to extend to the quantum domain the previous classical results. 
Our treatment, essentially algebraic in nature, is a re-elaboration and — we believe — a 
simplification of the technique originally proposed by R.G. Littlejohn in classical mechan- 
ics. It is based on a different choice of non-canonical variables adapting to classical as 
well as quantum mechanics. Depending essentially on the canonical structure the method 
applies indistinctly to the classical and the quantum theory. We nevertheless focus on the 
quantum problem. 

In order to better understand what is going on in the strong-filed regime of a quantum 
particle moving in an external magnetic filed it is better to first have in mind the main 
features of corresponding the classical problem [^]. Let us therefore to briefly consider 
a classical particle of mass m and charge e moving in a homogeneous magnetic filed of 
intensity B. As is well known the trajectory of the particle is represented by a spiral 
wrapping around a field line, as sketched in Fig.[^a: the particle performs a uniform 
circular motion of frequency ujb = eB/mc and radius rs = mc\v±\/eB {\v±\ is the norm 
of the normal component of the velocity) in the directions normal to the field, while the 
center of the orbit, called the guiding center, moves freely along a filed line. Keeping fixed 
the initial condition, the stronger the magnetic field the faster the rotation of the particle 
when compared with the drift along the filed direction and smaller the portion of space 
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Figure 1: Different time scale behaviours of a charged particle in a strong magnetic field: 
a) fast rotation of the particle and slow guiding center drift in a homogeneous field; bj in 
the inhomogeneous case the guiding center drifts away from the filed line very-slowly. 



explored by the particle around the filed line. This indicates, to the one side, the presence 
of different time scales in the dynamics of the system and gives, on the other hand, the 
reason why the motion in a very strong magnetic filed may be studied along the same 
lines as that in a weakly inhomogeneous one. Let us introduce now a small inhomogeneity 
in the field. The picture of the motion should not change substantially. The particle, in 
fact, keeps on rotating around its guiding center, the frequency and the radius weakly 
depending now on the position, and the guiding center still drifts along a field line. In 
this case, however, the guiding center do not remains exactly on a single field line, it 
starts drifting very-slowly in the directions normal to the filed. Three different time scale 
behaviours of the system may therefore be distinguished: the fast rotation of the particle 
around the guiding center, the slow drift of the guiding center along a magnetic field line 
and the very-slow drift of the guiding center in the direction normal to the field. The 
situation is sketched in Fig.|I]b. The stronger the magnetic field the cleaner the separation 
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among the three degrees of freedom. 

An outlook to the canonical structure of the homogeneous case makes immediately 
clear how the introduction of kinematical momenta and guiding center operators allows to 
separate the three degrees of freedom of the system. This is briefly reported in section 2 
where the relevant notations of the paper are also set up. After discussing the geometrical 
complications involved in the adiabatic motion of a charged particle in an inhomogeneous 
magnetic field, section 3, an appropriate set of non-canonical operators generalizing the 
one used in the discussion of the homogeneous problem is constructed in section 4. These 



are obtained as formal power series in the magnetic length Ib = yhc/eB which appears 
naturally as the adiabatic parameter of the theory. The Pauli Hamiltonian describing 
the motion of the particle is then rewritten in terms of the new adiabatic operators in 
sections 5 and 6, whereas the anholonomic effects appearing in the adiabatic separation 
of the fast and slow/ very-slow degrees of freedom are discussed in section 7. Our results 
are summarized in equations (0), (|4^) and (|4^). In the classical limit these reproduce 
correctly the classical theory. Section 8 contains our conclusions. 

2 Canonical structure of the guiding center motion 

Magnetic interactions appearing essentially as modifications of the canonical structure 
of a dynamical system it is worthwhile to start by briefly discussing this peculiarity in 
the elementary case of a quantum charged spinning particle in a homogeneous magnetic 
filed. This allows us to immediately focus on the heart of the problem establishing at 
the same time terminology and notations. We consider therefore a spin-1/2 particle of 
mass m, charge e and gyromagnetic factor g moving in space under the influence of the 
homogeneous filed B{x) = B z. As in the inhomogeneous case, to be discussed later on, 
the physical dimension of the magnetic field is reabsorbed in the scale factor B, the inverse 
square root of which, opportunely rescaled, will play the role of the adiabatic parameter of 
our theory. Introducing an arbitrary choice of the vector potential a for the dimensionless 
filed B{x)/B, rota = z, the motion of the particle is described by the Pauli Hamiltonian 



V = {dx,dy,dz) denoting the gradient operator and cr = {crx,cry,az) the matrix-valued 
vector constructed by means of the three Pauli matrices. As well known the solution of 
this simple problem was first obtained by Landau at the beginning of the thirties and leads 
naturally to replace the standard set of canonical operators Pi = —ihdi, x*, i = 1,2,3, 





(1) 
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by the gauge invariant kinematical momenta iXi = Pi — {eB/c)ai and the guiding center 
operators X* = + {c/ eB)e^^TTj. A very rapid computation yields the nonvanishing 
commutation relation among the new variables 



[vr2,vri] = -i 



MB 




. he 



(2) 



c 



indicating 712-^1, tts-X^ and X^-X"^ as couples of conjugates variables. Moreover, the scale 
dependence of the commutators (|]) allows to identify the three couple of operators as de- 
scribing respectively the fast, the slow and the very-slow degrees of freedom of the system 
(see eg. 0). In terms of the new gauge invariant operators Hamiltonian (||) rewrites in the 
very simple form Ti = (vti^ + 7r2^ + 7r3^)/2m + gheBa^/mc. The harmonic oscillator term 
(tti^ + 7r2^)/2m takes into account the rapid rotation of the particle around its guiding 
center while the free term n^'^ /2m the slow drift of the guiding center along the straight 
magnetic field lines. The very-slow variables X^ and being constant of motion, the 
guiding center do not move in the directions normal to the field. Let us stress that in the 
canonical formalism the spatial rotation of the particle around its guiding center is taken 
into account by the phase space trajectory of a couple of conjugate variables: the particle's 
velocity components in the directions normal to the field: tti and tt2- The presence of an 
external magnetic field produces therefore a rotation of the canonical structure, mixing 
up spatial coordinates and canonical momenta in new canonical operators adapting to 
the different time scale behaviours of the particle! In section 4 we will construct such 
"adapted operators" — as power series in the adiabatic parameter — for the motion of a 
quantum charged spinning particle in an arbitrary magnetic filed. This allows to extend 
to quantum mechanics the Hamiltonian approach to the classical guiding center motion 
developed by R.G. Littlejohn The case of a magnetic filed with constant direction has 
been previously considered in . Before to proceed some preparatory material is however 
necessary. 

First of all it is convenient to introduce dimensionless quantities by factorizing the 
energy scale hujB-, ojb = eB/mc, from the Hamiltonian. This leads to redefine kinematical 
momenta and guiding center operators as 



= Jhc/eB being the magnetic length. The relevant commutation relations may so be 



TTj = —ilsdi — Ib CLiix) 
X' = x' + lBS'^TTi. 



(3) 
(4) 
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recast in the compact and very convenient form 




fast 
slow 



very-slow 



(5) 



where the spin variables have also been considered. 

As a second and more serious task the geometrical structure responsible for the an- 
holonomic effects appearing in the adiabatic motion in a strong magnetic field has to be 
discussed. 

3 Magnetism and geometric- magnet ism 

The beautiful analysis of the adiabatic separation of fast and slow degrees of freedom in 
a quantum system proposed by M.V. Berry H. Kuratsuji & S. lida [§], J. Moody, A. 
Shapere & F. Wilczek p, R. Jackiw |T^ and others, has pointed out that in lowest order 
the reaction of the fast to the slow dynamics is through a geometry-induced gauge struc- 
ture resembling that of (electro-) magnetism. This phenomenon has been identified and 
found to be important in a variety of physical contexts |Tl[] and has been recently referred 
by M.V. Berry & J.M. Robbins as geometric- magnetism [|12|. A rather curious fact, first 
pointed out by R.G. Littlejohn in a series of beautiful papers on the canonical structure of 
classical guiding center motion is that, in some circumstances, magnetism itself may 
generate geometric-magnetism. The aim of the present section is that of discussing the 
geometry involved in such "magnetic-induced geometric-magnetism". 

For shake of clearness it is useful to begin by briefly recalling the geometrical character 
of the kinematical quantities characterizing the motion of a particle in space. This will 
led to a rather intuitive picture of the geometrical structure involved in the adiabatic 
motion of a charged spinning particle in a strong magnetic field, allowing, at the same 
time, to frame it in a general and rigorous context. As well known the state of a particle 
moving in space is completely characterized by its position x and its velocity v, i.e. by a 
point in the tangent bundle TB? of the three-dimensional Euclidean space B?. The flat 
parallel transport of B? makes it natural to parameterize every fiber T^R^ of the bundle 
by means of a fixed reference frame in i?^, that is, to identify the tangent space in every 
point X with the physical space itself. Such an identification is certainly very useful in 
most circumstances, but it is a convention after all. In principle we are free to choose 
arbitrarily the frame of T^R^ in every x, the parallel transport — and not the way in which 
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Figure 2: Framing the tangent bundle TR^ of the physical space: a) hj means of a single 
fixed frame in R^; bj by using local reference frames adapting to the magnetic field lines 
geometry. 

we describe it — being all that matters. This freedom of arbitrarily rotating the reference 
frame of the tangent space in every point x, a local SO (3) symmetry, plays a crucial role 
in what follows. To visualize the situation, therefore, we shall picture the Euclidean space 
as filled up with orthonormal reference frames. To start with, we can imagine all them 
as combed parallel to a single fixed frame {x, y, z} in R^ (see Fig.|^), but even in a flat 
geometry, this is not always the better choice. 

The magnetic line bundle 

As qualitatively sketched above, the motion of a charged-spinning particle in a strong 
magnetic filed is characterized by the separation of time scales in the three degrees of 
freedom, making the system amenable to a perturbative analysis. In the lowest order 
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approximation the particle performs a fast rotation in the plane normal to the field line 
at which its guiding center is located. This is taken into account by the two components 
normal to the field of the particle's velocity (to this order a couple of conjugate variables). 
Disregarding the slow drift of the guiding center along the filed line and the very-slow 
motion, therefore, the velocity of a particle which guiding center is located in x is ef- 
fectively constrained to the plane /x^ generated by the vectors normal to the field in x. 
In every point of the space the magnetic filed b{x) picks the complex line out of the 
tangent space TgR^, reducing the tangent bundle TE? to a complex line bundle, hereafter 
the magnetic line bundle Ai^ It is then natural to use the local S0{3) symmetry of 
the theory to adapt the parameterization of TR^ to the subbundle A4 by combing, say, 
the z direction of the frame of every T^-R'^ according to the direction of the field. We so 
smoothly introduce point dependent adapted reference frames {ei, 62, 63} in such a way 
that in every point ei(x), 62(0;) parameterize fj,g while 63(0;) is aligned with b{x) (see 
Fig2b). Such reference frames are commonly used in the discussion of geometrically non 
trivial physical problems such as in general relativity and are referred as anholonomic 
frames. It is worthwhile to note that fixing 63 according to the filed direction reduces the 
local S0{3) symmetry oiTR^ into the local 5*0(2) = U{1) symmetry of Ai. The vectors 
ei{x) and 62 (^) are in fact determined up to the rotation 

ei(f) ei(x) cosx(f) - e2(x) sinx(^) 
e2{x) ei(x) sin X (if) + 62(0;) cosx(a;) 

x{x) being a point dependent angle. This residual ambiguity will result in the gauge 
freedom of our theory. 

Magnetic line bundle geometry 

We may now wonder how the vectors lying in A4 are transported from point to point, 
that is, whether the geometry of the magnetic line bundle is trivial or not. To this task 
we proceed in two steps. Considering a vector w{x) = w^e^{x), z/ = 1,2, in /j,^, we first 
transport it from the point x to the infinitesimally closest point x + dx hj means of the 
Euclidean parallel transport of and, second, we project it onto the plane i-i^^dx- (i) the 
Euclidean parallel transport of w in x + dx may be immediately evaluated as 

w{x + dx) = w{x) — w'^ (cj, • dkei)dx^ e^, 

^This subbundle of TR^ may be identified with the plane bundle of B. Felsager & J.M. Leinaas jist . 
See also the related paper of F. GHozzi |p^ . 
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latin indices running over 1, 2, 3, greek indices over 1, 2 and where the sum over repeated 
indices is imphedQ. The three quantities^ ei • 5^62, ei ■ due^ and 62 • dkC^ characterize the 
flat parallel transport of B? in the anholonomic frame (it is in fact possible to make them 
vanishing by rotating the adapted frames {61,62,63} back to fixed directions in every 
point), (ii) the projection onto ^i^j^^^ yields then 

w{x + dx)\^ = w{x) — (61 ■ dke2)dx^e^ 

indicating that the parallel transport along the infinitesimal path connecting x to x + dx 
produces the vector w to be rotated by the infinitesimal angle da = (61 ■ dke.2)dx^. When 
parallel transported along a finite closed path 7 the vector will therefore return to the 



starting point rotated by the angle |13 



= j (ei ■ dke2)dx^-, 



a 

this quantity being in general different from zero, the geometry of the magnetic line bundle 
results not flat. The operation of locally projecting onto the plane /x reduces the trivial 
5*0(3) local symmetry of the theory to a non-trivial 5*0(2) = f/(l) local symmetry! This 
local structure is described by a magnetic-like U{1) gauge theory. The parallel transport 
of the magnetic line bundle TW, results in fact completely characterized by the vector 

Ak = ei- (9fc62, (7) 

the connection one-form oi AA.. A. appears in the theory as a geometry-induced vector 
potential (to not be confused with the vector potential a representing the real magnetic 
field b). A point dependent redefinition of the local basis {ei{x), 62(x)} plays in fact the 
same role of a gauge transformation, the rotation (|^) producing the vector (|^) to transform 
according to Ak ^ Ak + dkX- The associate geometry-induced magnetic filed Bk = 
£kmnJ3mn, = dmAn — dnAm the cuTvatuTe two-forw. of A4, may also be considered. It 
is obviously a gauge invariant quantity and, being the rotor of a vector field, satisfies the 
Bianchi identity divS = 0. 

While the geometry-induced vector potential A. completely characterizes the intrinsic 
geometry of the magnetic line bundle Ai, the other two quantities 

hk = 61 ■ 5^63, 
hk = 62 ■ dkG^, 



^This notation will be employed throughout the rest of this paper. 

^The vectors ei, 62 and 63 being orthonormal in every point x, • Bj = Sij, these are the only 
independent quantities. 
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describing the flat parallel transport of in the anholonomic frame {61,62,63} may 
be seen a taking into account its extrinsic geometry. Since the curvature of the tangent 
bundle TR^ is zero the three quantities A., l\ and I2 are obviously not independent, being 
related by the equivalent of the Gauss, Codazzi-Mainardi and Ricci equations. The ladder, 
as an example, allows to re-express the geometry-induced gauge field 13 entirely in terms 
of Zi and I2 as 

B = hAl2, (9) 

^ki = {hkhi — hihk)l'^i ^ indicating the external product of R^ [|13|. With respect to 
the point dependent rotation (H) Zi and I2 transform as vectors (Zi — Zi cosx — I2 sinx, 
I2 h sin X + ^2 cos x) making the gauge invariance of B manifest. 



Magnetic filed lines geometry 

Thought the geometry of a magnetic filed is completely characterized by two independent 
function (e.g. the two independent components of the real magnetic field 6, or of the 
geometry-induced magnetic field B, etc. ) it may be useful to look at the problem from 
different points of view. We may wonder, as an example, how the intrinsic/extrinsic 
geometry of the line bundle is related to the geometry of magnetic filed lines. To this 
task we start by observing that the projection along the field direction of the two vectors 
Zi, I2 may be identified with the two second fundamental forms of the embedding of the 



magnetic field lines in R |jT^ . In every point of the space the curvature k of the magnetic 



filed line going through that point may so be expressed as 



k = ^ies-hy + ie,-l2r. (10) 

In a similar way the projection along the filed direction of the geometry-induced vector 
potential A have to be identified with the normal fundamental form of the embedding of 
the field lines in R^ (i.e. with the connection form induced by the Euclidean geometry onto 
the normal bundle of every filed line) JTSf. Up to the gradient of an arbitrary function, 
representing again the freedom of arbitrarily rotating the reference frame in the normal 
planes, in every point of the space the torsion r of the magnetic field line going through 
that point may be written as 

r = 63 ■ A (11) 

Curvature and torsion completely characterize the geometry of every single magnetic filed 
line and contain, in principle, all the informations relative to the geometry of our problem. 
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On the other hand we may also wonder about the global properties of the foliation of 
in terms of field lines. Of particular relevance for the adiabatic motion of a particle 
in an external magnetic field is the possibility of foliating space by means of surfaces 
everywhere orthogonal to the field lines. By virtue of Frobenius theorem this is controlled 
by the vanishing of the scalar = 63 • rotes. In terms of the magnetic line bundle 
geometry 

T = e,l2-e2h. (12) 

The magnetic filed lines torsion r and the Frobenius invariant JF play a crucial role in 
the description of the anholonomic effects appearing in the adiabatic motion of a charged 
particle in a strong magnetic field. 

4 Adiabatic quantum variables 

We are now ready for the construction of a set of adiabatic operators adapting to the 
different time scale behaviours of a quantum particle in a strong, but otherwise arbitrary, 
magnetic field. Let us consider therefore a spin-1/2 particle of mass m, charge e and 
gyromagnetic factor g moving in space under the influence of the inhomogeneous magnetic 
field B{x) = Bb{x), the physical dimension of the filed being again reabsorbed in the 
scale factor B. Denoting by a an arbitrary choice of the vector potential, rota = b, the 
dynamics of the system is described by the Pauli Hamiltonian 

H/hujB = ^-n-iTTi + gbi{x)ai, (13) 

where the kinematical momenta tTj = —ilBdi — ai{x)/lB have been introduced. The 
inhomogeneity of the magnetic filed makes Hamiltonian (|T3|) to depend on the position 
operators x, explicitly through spin term g bi{x)ai and implicitly through the commutation 
relations of the tTjS. In spite of the simple quadratic dependence of ( [TB| ) on the kinematical 
momenta, tti and 7^2 are in fact no longer conjugate variables and neither commute with 
TT^: the set of operators {TTi,x'^;i = 1,2,3} fulfil the commutation relations 

['7Ti,TTj] = ibij{x), [ni,x^] = -UbSI, [x\x^]=0, (14) 

bij{x) = eijkbk{x) denoting the skew-symmetric two-form associated to the field. In the 
lowest approximation we nevertheless expect the relevant degree of freedom of the system 
to be taken into account by the two components of the particle's velocity normal to 
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the filed. Considering the position operators x*s as adiabatic parameters driving the 
fast motion of the system we expect therefore the rapid rotation of the particle around its 
guiding center to be separated from the slow and very-slow motion by simply referring the 
kinematical momenta to the adapted frames introduced in the previous section. For the 
shake of concreteness we shall indicate by Ri'{x) the point dependent rotation bringing 
the fixed frame {x,y,z} into the adapted frame {ei(s), 62(0;), 63(0;)}. This allows to 
decompose the field b{x) in terms of its norm h = y/b ■ b and its direction b/b = Ri' Zj as 
hi{x) = b{x)Ri-' {x)zj. Once the rotation has been performed the kinematical momentum 
along the field direction decouples, up to higher order terms in the adiabatic parameter Is, 
from the other two components. The commutator of these, on the other hand, results to 
be proportional to b{x). Stated in a different way, in the adapted frame the particle sees 
an effective magnetic filed of constant direction and intensity b{x). To make the velocity 
components normal to the filed in a couple of conjugate operators it is now sufficient to 
rescale them by the point dependent factor b~^^^{x) (see 0). We shall indicate by Di^(x) 
the point dependent dilatation Di^ = diag(6^/^, 6^/^, 1) rescaling the first and second 
components of a vector by 6^/^ and letting the third one unchanged. 

In order to construct operators adapting to the fast time scale behaviour of the system 
two point dependent operations have therefore to be performed: (i) a rotation Ri' {x) to 
the local adapted frame and (ii) a dilatation Di^{x) rescaling the normal components of 
the kinematical momenta. The particle coordinates being not external parameters but 
dynamical variables of the problem these operations will produce higher order corrections 
in the various commutators. We shall therefore proceed order by order in the adiabatic 
parameter Ib by constructing sets of adiabatic operators fulfilling the desired commutation 
relation up to a given order in Z^: at the n-th order we shall look for a set of operators 
{nj"\ ^{n)'i ^ = 1) 2, 3} fulfilling the conditions 

• ria"^ are a couple of conjugate operators up to terms of order Ib"', 

• 113"^ -^(n)' -^fn) commute with n^"'^ Ila"^ up to terms of order Ib", 

• in the limit of a homogeneous filed, b{x) z, the adiabatic kinematical momenta 
H^i^^s and guiding center operators X^^-jS should reduce to the expressions and 
(I) respectively. 

Our present task being that of separating the fast degree of freedom from the slow and 
very-slow motion, we do not insist for the moment X^^-^ — U^^^ and X^^-^ — X'^^-^ to be 
couple of conjugate operators as in the homogeneous case. 
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For computational proposes it is very convenient to use a compact notation which 
does not distinguish among the physical dislike directions along and normal to the field. 
This probably obscures a while the physical contents of the various expression but greatly 
simplifies formal manipulations. When necessary we will therefore expand the notation in 
order to shad light on the physics involved. For the moment we proceed in the opposite 
direction by introducing the point dependent matrix 

(3,Kx)^Df{x)Rf{x) (15) 

representing the successive application of the two operations necessary to construct the 
adapted kinematical momenta in the lowest order. This allows to rewrite the skew- 
symmetric two- form hij{x) in terms of Ski — Skiz (representing a homogeneous filed directed 
along z) 

h,{x)=^f{x)(5f{x)eu. (16) 

The matrix jSi-' and this representation of the filed result to be very useful in the con- 
struction of the adiabatic quantum variables. 



Zero-order operators 



In order to construct the zero-order operators fulfilling the desired conditions up to terms 
of order it is sufficient to operate the rotation and the dilatation discussed above 



n, 



(0) 



(17) 



the matrix f3i^ being evaluated in X{q) = x. The anticommutator { , } is obviously intro- 
duced in order to make the n-'^'*s Hermitian. A rapid computation confirms our deductions 
yielding the commutation relations fulfilled by the zero-order adiabatic operators as 



-i Ib 
0, 



• hkl 



where F^/ 



(18) 



n^'^'' and n2°'' are 



{dkl3i ^)I3~1-' and all the functions are evaluated in X(o) 
conjugate operators up to 0{Ib)- The commutators depend on the derivative of magnetic 
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filed through the vector-valued matrix 



2 b 



/ 



(19) 



allowing to clearly distinguish the effects produced by a variation of norm of the magnetic 
filed from that produced by a change of direction. The ladder are entirely geometrical in 
character being taken into account by the magnetic line bundle connection form A. and 
by the two extrinsic vectors l\ and ^2- 



First-order operators 



Whereas the construction of the zero-order operators is in some way suggested by the 
physics of the problem, a more technical effort is required for higher order terms. The form 
of the first-oder guiding center operators is nevertheless suggested by the corresponding 
homogeneous expression (H), 



X, 



(1) 



X, 



(0) 



{/?.\nr} 



(20) 



the matrix j3k being again evaluated in X(o). We immediately obtain the new commuta- 
tion relation 



^(l)'^(l) 



2 



hkl 



-ilB^e'^'Pk'Pi' + Oils'), 



(21) 



indicating the 0{Ib^) decoupling of the adiabatic guiding center operators from Yi^i^ and 
n2°'*. All the functions are now evaluated in X(i). Though our analysis will be carried out 
up to 0{Ib^)-i we also wrote the the first nonvanishing contribution to the commutators 
among the X(\)S, which is of order Ib^- Even if unimportant in the present calculation, 
this allows us to visualize the very-slow time scale of the system. 

The construction of the first-order kinematical momenta is performed by looking for 
order Is counterterms to be added to the n|°'''s. These should be homogeneous second 
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order polynomial in the nj-°'''s with coefficients depending on A rather tedious 

computation produces 



n 



(1) _ TfiO) I 7 Mmn 

- iij + Lb c^j 



(22) 



where d. 



klmn 



24-"'- + Spi)^"''" + l^f^'^'^i^j + Spi)ehge'^'^'' and all the functions are 
evaluated in When expanded these expressions do not look so complicated as a ffist 
sight. We nevertheless insist in keeping this notation which greatly simplifies the following 
manipulations. The commutation relations among the first-order adiabatic variables are 
obtained as 



^(i)'^(i) 




-tlB^e'^'Pk'Pi' + OilB') 



(23) 



It is very interesting to observe that a monopole singularity, that is a point of nonvanishing 
divergence, represents an obstruction in the construction of the adiabatic operators. Being 
concerned with real magnetic filed we nevertheless assume div 6 = and carry on in our 
adiabatic analysis. n[^'' and 112^^ are then conjugate operators commuting with all the 
remaining variables up to terms of order /^^ and the fast degree of freedom decouples 
from the slow and very-slow motion up to terms of this order. 



A non-canonical set of operators 

At least in principle it is possible to repeat this construction an arbitrary number of times 
getting, as power series in Ib, a set of adiabatic non-canonical operators {Ili,X'^;i = 
1, 2, 3} fulfilling the commutation relations 

[n„ n,] = t [n„ X^] = -i h 6f Rf, [X\ X^] = -i e^' Rtnr, (24) 

all the functions being now evaluated in X. These formal series are in general — and have 



to be |T6[ — not convergent, representing anyway a very useful tool in the discussion of the 
adiabatic behaviour of the system. The description of the problem to a given order n in 
the adiabatic parameter Ib requires the knowledge of the first n + 1 terms of the adiabatic 
series, so that up to terms of order Ib^ we may identify the IljS and X*s with the Ii^l^?> 
and -^(i)S respectively. An outlook to the commutation relation (24) allows to clearly 
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identify the dependence of the canonical structure on the variation of norm and direction 
of the magnetic filed. Whereas a suitable redefinition of reference frames in TR^ allows 
to separate the fast degree of freedom from the others, the very-slow variables are made 
into a couple of non-conjugate operators by an inhomogeneous intensity while a variation 
of the filed direction even produces the mixing of very-slow and slow variables. The 
description of these by means of couples of conjugate operators requires the introduction 



of curvilinear coordinates in space [0 , the so called Euler potentials |18[ . We do not insist 
further on this point for the moment observing that under the action of 111,112 and lis, 
X^, X^, the Hilbert space of the system separates in the direct sum of two subspaces 
describing respectively the rapid rotation of the particle and the guiding center motion. 



5 Expanding the Hamiltonian 

The adiabatic operators 11 and X constructed in the previous section have been introduced 
in such a way to embody the expected features of the motion of a quantum charged particle 
in a weakly- inhomogeneous magnetic field. Their main advantage lies, in fact, in the very 
suitable form assumed by the Pauli Hamiltonian when rewritten in terms of them. To 
this task we have first to invert the power series expressing Ilj and X* in terms of the 
operators vTjS and x*s and, second, to replace these in (|1^). This yields the Hamiltonian 
as a power series in the magnetic length Ib, 

n = n^''^ + iBn^^^ + iB^n^^^ + ... (25) 

allowing the adiabatic separation of the fast degree of freedom from the slow/very-slow 
motion and the evaluation of approximate expressions of the spectrum and of the wave 
functions of the system. In order to get the vTjS and x*s in terms of the HjS and X*s we first 
recall that X^ = X^-^-, + 0{Ib^)- By rewriting X^^-, in terms of the Tlf \ and X^q-, = x^s, 

(pop, Yif'^ in terms of the vr^s and x*s, (p!7|), and by solving with respect to x*, we then 
obtain as a function of the vTjS and the X*s, x* = x*(7f, X). This allows to rewrite H 



(0) 
i 

as a function of the vTjS and X*s. Recalling finally that Hj = Tl^P + 0{Ib^) and using ( ]22|) 
we immediately get Hj in terms of the tt^s and X*s, Hj = Hj(7f, X). The inversion of this 
relation, order by order in Z^, allows to get tTj and x* in terms of the adiabatic operators. 
The computation gives 

= + ^^T''[P-'iP^^on.{^k,Tll]] + 0{Ib% (26) 

X* = X*-/b£'='/5,*H, + 0(/b'), (27) 
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where cfj'"" = ^S^SjE"^'' — 20^]™"". As a useful check the commutation relations (0) may 



be reobtained by means of the (^41). 

The substitution of (^), (0) in the Pauli Hamiltonian ( ]IB| ) yields immediately the 
first two terms of the adiabatic expansion (p5|), 



H(°V^^b = ^p-'^p-'lUkUi+gkai (28) 
n'^'^/huB = p-'[p-'rCprf^^^o^{^^^'}^i}-9e''f3kidhk)a.Ui (29) 

where the notation 111^11^11; | = flfcllgn; + n^ngllfc has been introduced. In order to get 
some more physical insight in this expressions we now abandon our compact notation 
in favour of a more transparent one. By recalling the definition (|T5[) of j3i-'{x), (p!9| ) 
of and r^j^{x) and the explicit expression of the inhomogeneous dilatation Di-'{x) = 
diag(6^/^(x), 6^/^(x), 1), we rewrite everything in terms of the magnetic field and of other 
quantities capable of a direct physical interpretation. The full expansion of the zero order 
Hamiltonian (^) gives 

H^^^/hooB =^Us^ + b[J + g{es- ct)] , (30) 

where J represents the harmonic oscillator Hamiltonian constructed by means of the 
canonical variables Hi and H2, J = (Hi^ + H2^)/2, and the norm of the magnetic field 
b{X) is evaluated in the adiabatic guiding center operators X. We observe that while the 
H1-H2 degree of freedom decouples to this order from the slow and very-slow variables 
the spin does not. The separation, up to higher order terms, of the fast motion (rotation 
+ spin) requires in fact a subsidiary zero order transformation which we will perform 
in the next section. For the moment let us observe that, up to the spin term, the zero 
order Hamiltonian ( pOf ) precisely embodies the expected behaviour of the system: the 
canonical couple of operators H1-H2 takes into account the fast rotation of the particle 
around its guiding center, while the non-canonical variables U^-X^ describe the slow 
motion along the magnetic field lines by means of an effective "kinetic energy -|- potential 
energy" Hamiltonian. The norm of the magnetic filed b{X) plays the role of the effective 
potential. As long as 0{Ib^) terms are ignored the very-slow dynamical variables X^-X"^ 
appear only as adiabatic parameters driving the slow motion, whereas a more accurate 
analysis indicates them as taking into account the very-slow drift in the directions normal 
to the field M. 
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More complicated appears the full expression of the first order Hamiltonian (^) . The 
replacement of (5i'{x) and V^^^i^x) by means of (|15D and ([T9|) yields in fact the expression 



-b'/\e,.A)J, 



+ 



-(ei ■ Za - 62 ■ Zi) - (eg ■ A) 



JUs + -(ei ■ Z2 + 62 ■ Zi) (Hi^ - U2^) 



4 



■-(ei • Zi - 62 • Z2) {Hi, 02} Hs + b-'/^ [(ea ■ Z2) Hi - (eg ■ Zi) Ha 



-gb'/^e^-'Ue^ ■ Zi)(ei ■ a) + (e^ ■ Z2)(e2 ■ a] 



(31) 



indicating the first order coupling among the various operators. The notation = 
|5"^nQn^n/3 has been introduced and all the functions are evaluated in X. As expected 
from dimensional considerations Ti.^^^ depends only on the first order derivatives of the 
filed. It is nevertheless worthwhile to stress that the gradient of the magnetic-field-norm, 
grad b = V6, appears only in the first term of the right hand side of this expression, all the 
remaining terms depending only on the quantities A, li and Z2 completely characterizing 
the intrinsic/extrinsic geometry of the magnetic line bundle Ai. To a large amount, 
therefore, the complication of this expression is produced by the variation of direction of 
the magnetic field, that is, by the nontrivial geometry of It is not yet time to comment 
on the structure of Ti.^^^ First of all, it is in fact necessary to operate a suitable unitary 
transformation separating the zero order fast motion from the other degrees of freedom, 
that is diagonalizing the the spin term 63 ■ cr. This will produce a modification of the 
first order term of the adiabatic expansion. Secondly, it is possible to drastically simplify 
the form of Ti*-^-* by operating a suitable first order unitary transformation. The strategy 
is nothing else than the quantum equivalent of the so called averaging transformation of 
classical mechanics and results of great help in shading light on the physical content of 



6 Quantum averaging transformations 

A well known strategy in dealing with the adiabatic separation of fast and slow variables in 
classical mechanics consists in performing a series of successive canonical transformations 
(the averaging transformations) separating, order by order in some adiabatic parameter, 
the rapid oscillation of the system from its slow averaged motion. The analysis depending 
essentially on the canonical structure of the problem generalizes immediately to quantum 
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mechanics, the canonical transformations being replaced by suitable unitary transforma- 
tions. The full adiabatic expansion describing the motion of a spin degree of freedom 
adiabatically driven by external parameters has been obtained along these lines by M.V. 
Berry while R.G. Littlejohn and S. Weigert employed the method in discussing 



the first adiabatic corrections to the semiclassical motion of a neutral spinning particle in 
an inhomogeneous magnetic field. We shall consider therefore a set of unitary operators 

?7(") =exp{i/B"/:(")}, (32) 

n = 0,1, ... such that fast and slow/ very- slow degrees of freedom separate up to 0{Ib"'^^) 
in the Hamiltonian obtained by the successive application of U^^\ U^^\ ...,[/("). Whereas 
in classical mechanics it is natural to consider the averaging transformation as defining 
new canonical variables, in quantum mechanics it appears more convenient to keep the 
canonical operators fixed and transform the Hamiltonian. 

Zero-order transformation 

The zero order separation of the fast and slow/ very- slow motion requires the diagonal- 



ization of the spin term gh{X){e^[X) ■ cr) of Hamiltonian (30). Denoting by Pi{x) the 
infinitesimal generator of the rotation Ri'{x) bringing the fixed frame {x,y,z} into the 
adapted frame {ei{x) , e2{x) , e^{x)} , Ri^ = (e^)/ = 6j + pi^ + \pi^ Pk^ + the aim is 
achieved by choosing 

= -\e''''p^,{X)ak, (33) 

the matrix pij = pi^ being evaluated in the guiding center operators X. Because the 
commutation relations ( p^ the operator U^^^ commute with Hi, H2 and therefore with J, 
produces 0{Ib) terms when commuting with H3 and 0{Ib^) terms when commuting with 
functions of X. In evaluating the new Hamiltonian H' = U^^^HU^^^^ = 7^(°)' + /bH(^)' + ... 
up to terms of order we have therefore to worry only about the action of U^^^ on a 
and H3. A very rapid computation yields the transformation rule 

f/W(ei-o-)f/(°)^ = a, + 0(/B') (34) 

while the action of [/(°) on H3, [/(o^Hst/W^ = Hg + f/^^^pg, f/^^], may be easily evaluated 
by computing the commutator in the original set of operators tTjS and x*s and transforming 
back to adiabatic variables 

f/(o)n3f/(o)t ^ ng + hie, ■ h)a, - hie, ■ h)cx2 + hie-s ■ A)a, + 0{Ib'). (35) 
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Subjecting and 7^^^) to the zero order averaging transformation U^'^^ and by using 
(|3^ ) and (^Sj) we obtain the new adiabatic expansion 



(36) 



+ 



-(ei ■ Za - 62 ■ /i) J - (eg ■ A){J - 0-3) 



ns + 7(ei ■ /2 + 62 ■ /i) (Hi^ - Ha') 



1 



--(ei • /i - 62 • ^2) {Hi, 02} + r^/^ [(es • Z2) Hi - (eg ■ /i) 



+ 



[63 ■ 12)01 - (63 ■ Zi)(T2 



n3 - <^6V2^A^4(e^ . Zi)ai + (e^ ■ ^2)^2 



(37) 



All the functions are evaluated in X. The fast and slow/ very- slow motion are separated 
in this way in the zero order term of the adiabatic expansion but not in the first order 
term. 



First-order transformation 

The application of the first order averaging transformation U^^^ to 7i' produces the new 
Hamiltonian U" = U^'^'^H'U^^^^ = H^^^ + /^(^(i)' + i[C^^\n^^y]) + ... . It is then possible 
to simplify the first order term of the adiabatic expansion by choosing C^^'^ in such a way 
that its commutator with Ti,^'^^ cancels as much terms as possible of Ti'-^^ . The analysis 
of the commutation relation involved and a little thought indicates that it is possible to 
annihilate all but not the third term of (|37|) by choosing 



£(1) 



-6-3/2(e^. V6) + ga,U,^ +h-'/'6'^^ie,- A) J, 



-\h-\e, ■ ^2 + 62 ■ /i) {Hi, Hs} n3 - \h-\e, ■ h - ■ h) (Hi^ - Y[^^) 



(ea ■ h) n2 + (63 ■ Zi) Hi 



g b (63 ■ /2)cr2 + (63 ■ li)ai 



n. 



+ 



9 



9' 



[e^ ■ Zi)(2ai5^'^ - ga^e^^n - (e,, ■ Z2)(2a25^^ + gcr^e 



n„ 



(3J 



The commutators of the zero order Hamiltonian (^) with the various terms of C^^^ yields 
the terms of (|37D times the imaginary factor i, in such a way that they cancel in the 
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new adiabatic expansion. All the terms but not the third. It is in fact immediate to 
convince that no operator may be found in such a way that its commutator with (|36|) 
produces a term proportional to JII3 and cr^Us. The third therm of (0) may therefore not 
be removed from the adiabatic expansion representing a real first order coupling among 
fast and slow/very-slow motion and not a complication produced by a wrong choice of 
variables. Its relevance in the context of the classical guiding center motion has been 
first recognized by R.G. Littlejohn It is therefore not a surprise to re-find it in the 
discussion of the quantum guiding center dynamics. The quantum averaging method 
produces so the adiabatic expansion 

= (39) 



1 

-(ei ■ ^2 - 62 ■ li)J - (63 ■ A)(J - 0-3) 



n3, (40) 



We observe that whereas the zero order terms ( p6D depends only on the magnet ic-filed- 
norm b (other than on the commutation relations (p^)) the first order term (|iO|) is com- 
pletely characterized by the Frobenius invariant (0), and by the magnetic filed lines 
torsion (p!T|). 



7 Quantum guiding center dynamics and 
magnetic-induced geometric- magnetism 

The construction of a suitable set of non-canonical operators embodying the classically 
expected features of the motion of a charged particle in an inhomogeneous magnetic field 
and the quantum averaging method allow us to rewrite the Pauli Hamiltonian (p!3| ) in 
such a way that the fast degree of freedom — corresponding to the classical rotation of 
the particle around its guiding center — and the spin degree of freedom separate, up to 
terms of order Z^^, from the guiding center dynamics. The transformation to the adiabatic 
operators IljS, X*s, (^Dj) and (P^, and application of the zero and first order quantum 
averaging operators, (^) and (|38|), produces in fact the Hamiltonian 



n/hcUB = -U3^ + b{J + g(T3)-l 



B 



T{J-a,) 



2^. 



II3 + 0{Ib' 



(41) 



Disregarding terms of order higher than Ib the operators J, representing the magnetic 
moment of gyration of the particle, and are constant of motion of the system. Frozen 
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the particle in one of its J and as eigenstates Hamiltonian (^Tj) describes therefore the 
corresponding guiding centre dynamics. As long as 0{Ib^) are ignored and X"^ appear 
as non-dynamical external adiabatic parameters and only the Tl^-X^ degree of freedom, 
representing in the classical limit the drift of the particle along the magnetic field lines, 
is dynamically relevant. To this order, therefore, the quantum guiding center dynamics is 
described by a one degree of freedom Hamiltonian given by the sum of the kinetic energy 
113^/2 and of an effective potential proportional to h{X). lis being a slow variable, that 
is of the same magnitude of the first adiabatic correction, the order Is term [r(J — (T3) — 
jFJ/2]n3 may be identified as a magnetic-like interaction and reabsorbed in the zero order 
Hamiltonian as a gauge potential. The guiding center Hamiltonian rewrites in this way 
in the familiar form 

n/huB = ^ (Ha - Ib A{X)f + V{X) + 0{Ib^), (42) 

with 

A{X) = (J-a3)r(X)-^^(X), (43) 

V{X) = {J + gars)b{X). (44) 

As it might be expected from the general discussion of section 3 the magnetic filed line 
torsion t = e-^- A. appears as (a part of) a gauge potential in the effective slow dynamics, 
taking into account the anholonomy produced by the non trivial parallel transport of the 
magnetic line bundle A4. Maybe unexpected, at least form this point of view, is the 
contribution given by the Frobenius invariant. Let us in fact compare the guiding center 
motion of a charged particle along a magnetic filed line with the propagation of light in a 



coiled optical fiber |T9[ or to the motion of an electron constrained on a twisted line |^ 
In both cases — sharing the same mathematical background — the adiabatic dynamics is 
coupled with an effective vector potential proportional to the torsion of the optical fiber or 
of the twisted line. The analogue contribution appears in the guiding center motion, (J — 
(T3) t{X), but it is not the whole story. The particle being not homogeneously constrained 
in the neighborhood of a line, it results sensible to the variation of the geometry in the 
magnetic field lines surrounding the one on which the guiding center si located. If all 
the field lines would have the same geometry, the foliation of in terms of them would 
be trivial, the Frobenius invariant zero and the situation analogue the examples above. 
The geometry of this foliation being in general non trivial it yields a further contribution 
to the gauge potential A{X) proportional to the Frobenius invariant, JjF(X)/2. It is 



23 



obviously not possible, in the general case, to remove the gauge potential ( ^3] ) by means 
of a suitable choice of gauge. 

In order to make the identification of with a one — two, if we want to consider X^, 
X"^ as dynamical — degree of freedom Hamiltonian complete it is necessary to replace the 
IljS, X*s by a set of canonical operators. The task is achieved by introducing a Darboux 
coordinate frame x* = x*(a;), i = 1,2,3, bringing the closed skew symmetric two-form 



hij{x) in the canonical form. 



bki{x)-;--Tr- = ^ii (45) 



(x^(x), x^(a;) may be identified with a pair of Euler potentials for the magnetic filed 
b{x), Vx^ A Vx^ = b, while x^{x) with the arc length of the magnetic filed lines). Defining 
X* = x*(X), i = 1,2,3 we get the canonical commutation relations [n3,X*] = —UbSI and 
[X^X-^] = —Ub^s^^ allowing the identification of the operators describing the slow and 
the very-slow degrees of freedom. It is in principle possible to start from the beginning 
by introducing such curvilinear coordinates in and work out the problem by using 
a canonical set of operators |l^.0 Nevertheless, whereas the existence of a Darboux 



coordinate frame is always guaranteed by Darboux theorem, it is hardly ever possible 
to find it explicitly and to proceed to the construction of the X*s. For this reason — 
thought the IljS, X*s appear as the most natural variables for the problem — the explicit 
construction of a set of non canonical operators appears as a better strategy. 



8 Conclusions 

The main difficulty in addressing the separation of fast and slow degrees of freedom in 
the study of an adiabatic system consist generally in finding out a suitable set of variables 
adapting with sufficient accuracy to the different time scale behaviours of the system. 
Starting from the homogeneous case and the canonical commutation relations (^) we 
showed how the analysis of the canonical structure of a charged spinning particle moving 
in an external inhomogeneous magnetic field leads naturally to the construction — as power 
series in the magnetic length Ib — of a suitable set of non-canonical operators allowing to 
systematically take into account the coupling between spatial and spin degrees of freedom. 
The new variables fulfil the very compact commutation relations (|2^ ) clearly displaying 

^The introduction of Darboux coordinate would produce automatically the framing of TgR^ by means 
of an adapted frame {ei(a;), e2{x), 63(2;)}. Our method, on the other hand, consists in adapting the frame 
of TsR^ without introducing the curvilinear coordinates. 



24 



the dependence of the canonical structure on the norm and direction of the external 
magnetic field. In terms of the new operators the Pauli Hamiltonian rewrites as a power 
series in the adiabatic parameter Ib which may be brought in a particular simple form 
by operating suitable unitary transformations. In this way the fast degree of freedom of 
the system representing classically the rapid rotation of the particle around the guiding 
center and the spin are separated from the remaining degrees of freedom up to terms of 
order The resulting effective guiding center dynamics displays geometric- magnetism: 
the coupling with the geometry induced gauge potential depending on the magnetic 
filed lines torsion ([Tl| ) and on the Frobenius invariant (|12|) , and with the scalar potential 
P^), proportional to the magnetic field norm. This completely extend to the quantum 
domain the previous classical treatments of the problem showing that the anholonomy 
first studied by R.G. Littlejohn in the classical guiding center theory plays an equivalent 
role in the discussion of the quantum problem. It is a feature of the canonical structure 
of the system after all. The geometrical mechanism responsible for the appearance of 
induced gauge structures has also been analyzed in some detail and formalized in the 
geometry of the magnetic line bundle A4. 

In concluding we observe that our discussion gives in some sense the solution of only 
half of the problem. The guiding center dynamics is still characterized by the presence of 
a fast and a slow time scale {slow fast, very-slow —* slow) and is therefore amenable 
to a treatment by means of adiabatic techniques. Nevertheless, the remaining problem is 
not of a so deep geometrical nature as the original one and is probably not capable of a 
treatment in general terms. 

Acknoledgments 

It is a genuine pleasure to thank M.V. Berry for hospitality and a very stimulating "few 
days discussion about fast and slow'^ in Bristol. I am also indebted with J.H. Hannay, E. 
Onofri and J.M. Robbins for very helpful discussions on related topics. 



25 



References 

[1] R.G. Littlejohn, J. Math. Phys. 12, 2445 (1979); Phys. Fluids. 29, 1730 (1981); 

J. Plasma Phys. 29, 111 (1983); in Contemporary Mathematics Vol.28, 151, J.E. 
Marsden ed. , (American Mathematical Society, 1984); Phys. Rev. ASS, 6034 (1988) 

[2] M.V. Berry, in Fundamental Aspects of Quantum Theory, 267, V. Gorini and A. 
Frigerio eds. (Plenum, New York, 1986) 

[3] Y. Aharonov and A. Stern, Phys. Rev. Lett. 69, 3593 (1992) 

[4] R.G. Littlejohn and S. Weigert, Phys. Rev. A48, 924 (1993) 

[5] T.G. Northrop, The Adiabatic Motion of Charged Particles (Interscience, New York, 
1963) 

[6] R Maraner, J. Phys. A: Math. Gen. 29, 2199 (1996) 

[7] M.V. Berry, Proc. Roy. Soc. London A392, 45 (1984); in Geometric Phases in 
Physics, A. Shapere and F. Wilczek eds. , (World Scientific, Singapore, 1989) 

[8] H. Kuratsuji and S. lida. Prog. Theo. Phys. 74, 439 (1985) 

[9] J. Moody, A. Shapere and F. Wilczek, Phys. Rev. Lett. 56, 893 (1986) 

[10] R. Jackiw, Comm. At. Mol. Phys. 20, 71 (1988) 

[11] A. Shapere and F. Wilczek, Geometric Phases in Physics (World Scientific, Singa- 
pore, 1989) 

[12] M.V. Berry and J.M. Robbins, Proc. Roy. Soc. London A442, 641 (1993); Proc. Roy. 
Soc. London A442, 659 (1993) 

[13] B. Felsager and J.M. Leinaas, Nuc. Phys. B166, 162 (1980) 

[14] F. Gliozzi, Nuc. Phys. B141, 379 (1978) 

[15] M. Spivak, A Comprehensive Introduction to Differential Geometry, Vol.4, (Publish 
or Perish, Houston, 1975) 

[16] M.V. Berry, Proc. Roy. Soc. London A414, 31 (1987) 



26 



[17] C.S. Gardner, Phys. Rev. 115, 791 (1959) 
[18] D.P. Stern, Am. J. Phys. 38, 494 (1970) 

[19] R.Y. Chiao and Y.S. Wu, Phys. Rev. Lett. 57, 933 (1986); A. Tomita and R.Y. Chiao, 
Phys. Rev. Lett. 57, 937 (1986) 

[20] S. Takagi and T. Tanzawa, Prog. Theoret. Phys. 87, 561 (1992) 



27 



